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Abstract 

The simplest type III seesaw model as originally proposed introduces one lepton triplet. It 
thus contains four active neutrinos, two massive and two massless at tree level. We determine 
the radiative masses that the latter receive first at two loops. The masses are generally so tiny 
that they are definitely excluded by the oscillation data, if the heavy leptons are not very heavy, 
say, within the reach of LHC. To accommodate the data on masses, the seesaw scale must be as 
large as the scale of grand unification. This indicates that the most economical type in model 
would entail no new physics at low energies beyond the tiny neutrino masses. 
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1 Introduction 



The standard model (SM) of electroweak interactions when viewed as an effective field theory at 
low energies, has a unique dimension five operator that can generate Majorana neutrino masses 

And the operator has only three possible realizations at tree level [|2l. These correspond to 
the celebrated three types of seesaw models |l3l|4l[5]|. While the type I model introduces sterile 
neutrinos as the minimal option to operate the seesaw, the other two prescribe particles that 
participate electroweak interactions. If the seesaw scale is not too high, richer phenomena are 
expected in the last two types of models. There have been extensive investigations on the type I 
and II seesaw models, but the interest in type III has been catalyzed recently by the advent of the 
LHC at CERN [6], where the assumed triplet leptons could be directly produced through gauge 
interactions if they are not too heavy [l7l[8l|9lllOl. Various other phenomenological aspects of 
the model have also been explored, including possible modifications to leptogenesis [[TTlfTllfTSl 
[MlfTOl, low energy effects in lepton flavor changing processes ifTSl and anomalous magnetic 
moments of charged leptons lfT6l [TTll . renormalization group running of neutrino parameters 
ifTSl . and the potential role as dark matter [[T9ll . to mention a few. 

For a seesaw model like type III to be relevant at relatively low energies, it must be capable 
of incorporating the data from oscillation experiments and other constraints with a not too high 
seesaw scale. We are thus motivated to start with the simplest type III seesaw as was originally 
proposed |l5l. It extends SM by one triplet of leptons, resulting in two massive and two massless 
neutrinos at tree level, plus a pair of heavy charged leptons. It also serves as an approximation 
to more general structures that contain additional sequentially heavier triplets of leptons. The 
massless neutrinos not being protected by any symmetry should receive radiative masses, which 
will be determined in this work. It would be interesting to ask whether it is possible in this 
minimal model to get a radiative mass at a desired level with a seesaw scale accessible at LHC. 

The idea of generating a one-loop radiative mass for neutrinos was originally suggested in 
Ref EOl . and extended to two loops in [[2Tll22| . It offers a nice way to induce hierarchical and 
tiny neutrino masses. There is a vast literature that extends the idea in various aspects (see as 
examples, Il23ll2l[25l|26ll27l) and calculates radiative masses in different models [l28l[29l[30]l . 
We would not attempt to review the topic but reemphasizing the point that for a mechanism of 
radiative mass generation to be testable at colliders |I3T1 Ell [33l the relevant heavy mass scale 
cannot be too high. 

The paper is organized as follows. We describe in some detail the minimal model in the next 
section to set up our notations. The exact constraints on the lepton masses and diagonalization 
matrices are highlighted. They will be extensively utilized in our analytic evaluation of radiative 
mass. Also listed are the Yukawa couplings of leptons that may be useful in other applications. 
The radiative mass is then calculated in section [3] in a manner that facilitates later numerical 
analysis, and the final answer is given in terms of some loop integrals. These integrals are 
defined in Appendix A, and their leading terms in the heavy mass limit are given. For numerical 
analysis in sectionHl we first demonstrate the order of magnitude of radiative mass for a heavy 
mass scale that would be accessible at colliders. Then we consider the heavy mass limit trying 
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to accommodate neutrino masses derived from oscillation experiments. We conclude in the last 
section where the main points of the work are recapitulated. 



2 Type III seesaw model 

We describe systematically in this section the type III seesaw model proposed in Ref. [|5]|. While 
the exposed relations among the lepton mixing matrix and the lepton masses will be employed 
in the next section to evaluate the radiative neutrino masses, the displayed interactions may also 
be useful in other applications. 



2.1 Yukawa couplings and lepton mass matrices 

The model introduces a lepton multiplet, E, that is a triplet of SU {2)i but carries no hypercharge, 
on top of the fields present in SM. We shall restrict ourselves to the leptonic sector of the model. 
The lepton fields are 

\ / J_yO 



We have assumed without loss of generality that E is right-handed (RH). The Yukawa couplings 
plus the bare mass for E are 

-^Yuk = ^tr(MzE^E^ + M|E|E/,) 

+ {Ky^fR^ + ^^TrjIfi)! + {FlyjXR^ + ^^^Ry^Ft) , (2) 

where 4> is the scalar doublet with 4> = io^^* . y<i> and are respectively 3x3 and 3x1 
complex Yukawa coupling matrices. The superscript C denotes the charge conjugation, i//- = 
<^yO)^* with ^ = i'f''f'. Our notation is such that i//^ = {^fL)'^ ■ It is not necessary to include a 

— E^ coupling since ^f^X^ = X¥- Note that we can choose Mz, which is the seesaw scale 
in the model, to be real positive as any phase of it may be absorbed into yi,. 

After 4> develops a vacuum expectation value, v, the lepton mass terms become 



= ^Ms(EOEOC + E^E+c + E+E^c + h.c. 



+^ (^fLy<t>fR + -^^y^^R + fLyi-^R + h-c-) (3) 

Since E^ carry electric charge, they cannot be Majorana particles. Instead, their equal bare 
mass suggests the combination to a Dirac field, 

^ = E^+E+c (4) 
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with L^'^ = '^'/'(E^)*. It is then impossible to assign a lepton number to ^ without explicitly 
breaking gauge symmetry. The lepton mass terms are summarized as 

1 r 

-^m = -NLmNN^ + ELmEER+h.c. (5) 
where the neutral and charged lepton fields and their mass matrices are 



2.2 Gauge couplings of leptons 

The kinetic term for the triplet field is 

^4 = trL^ipLR, (8) 

where the co variant derivative is 

D^Zr = d^LR-ig2^[A"^a",i:R] (9) 

with and g2 being the SU {2)l gauge fields and coupling. The kinetic term can be ex- 
pressed in terms of the fields defined in eq Q. In so doing, the following relations are useful, 
yfCyt^^c ^ _^y^i^f^ ^ = W ~ ^ ^lOCY^ w) ^ where the total derivative may be dropped 
from Lagrangian. 

Including the standard kinetic terms for the SM fields Fl and fR, the complete kinetic terms 
for leptons are 

^kin = Ui^N + Ei^E + ^ (jw^W+ + ) + —J^Z^+eJ^A^, (10) 

Z V ^ / cw 

where , and axe the weak and electromagnetic fields coupled to the currents 

4^ = Nf{wLPL + WRPR)E 

Jz = Nf^i^PtN + E'f{z^LPL + ziPR)E 

JL = -Ey^E (11) 

and 7^y^ = {J^^)\ with the coupling matrices being 







/ (-5+4)13 \ E_f 



V 2 ■ -^M/.^. , --w^. (12) 



We have used the conventional notations cw = cos 6w, sw = sinOw, with 6w being the weak 
angle. 
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2.3 Diagonalization of lepton mass matrices 



Noting that the upper-left 3x3 block of m^f is zero, we can make standardized as follows. 
A unitary transformation in family space, Fi '^'^Fi, only modifies the Yukawa couplings, 
yz ^yi. and 3;$ ^y<i>. One can choose ^ to rotate the column vector to its third 
component, so that 

/O2 \ 
mN=\ ivrs (13) 

\ ivrs Ms / 

where rj: is real positive. There are thus two massless neutrinos (named 1 and 2) at tree level. 
They will generally get a radiative mass as their masslessness is not protected by any symmetry. 
The other two neutrinos (3 and 4) get the masses 

1 



/M2 + (vrE)2^M£ 
The mass eigenstate fields of neutrinos are therefore 



(14) 



where 



Vl = UI^Nl, Vr = v): = uINr (15) 



Un= I ice se | (16) 

-is0 C0 



with Ce = COS0, S0 = s'mO and tan0 = a/'«3/w4. 

The mass matrix of the charged leptons is diagonalized by bi-unitary transformations, 

El,r = Ui^dt^R, uIiueUr = diag(me, m^,mT,m;^;) (17) 

Here V4 and x are the new neutral and charged leptons beyond SM. They must be very heavy 
to evade the experimental detection so far. The tiny (small) mass of the observed neutrinos 
(charged leptons) then implies that, to very good precision, we have approximately 

m^^My, ^ — , (18) 
which will be employed in later numerical analysis. 



2.4 Summary of lepton interactions 

We can now express the interactions of leptons in terms of their mass eigenstate fields, V/ (/ = 
1,2,3,4) and (oc = e.jl, T^X)- The currents in eq. (fTTj) become 

JL = -ifi (19) 
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where 

= uIwlUl, Wr = UIwrUr 

= ulzfUL, = uIz^^Ur (20) 

Note that there is a degree of freedom in presenting the neutral current of Majorana neutrinos. 
Using V = Vl + vf = v*^ and y/^ Pl,rX'^ = -XY^Pr.lW^ we can write 

vy^^lPLV = \vy^ {3f^PL~^^^PR) V (21) 
Since f/^v and w, z are known, the following explicit results are useful: 
I2 

Wl = I ice -iV2s0 I Ul, Wr = V2 

S0 \flCQ 
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^l = A 4 ^CQSQ (22) 



-iCQSe 4 

One observes from the above that the right-handed charged current involves only the massive 
neutrinos V3^4 while the flavor changing neutral currents occur for both charged leptons and 
(massive) neutrinos. 

For completeness, we present some additional results that may be useful in other applica- 
tions of the model. First of all, one can construct the coupling matrices in the neutral currents 
in terms of those in the charged currents: 

^/ = I4 - = 4I4 - ^^z^^^L, = - \lfl'WR (23) 
The Yukawa couplings of the would-be Goldstone bosons G^'^ are 

= +-|^G+v[mv(#iPL + W)-(^LPi? + WL)m£]£ + h.c. 
\/2my^ 

G^v [my TIPl - ^l^vPr] V 

(24) 



cwmz 



m ^0 



mei^lPL + ^^Pr) - {S^tPR + T;,PL)m, 
Cwmz I J 

where my and are the diagonal mass matrices of the neutrinos and charged leptons. The 

above simple structure is dictated by the nature of G^'^ although the intermediate results in a 

direct derivation from .Sfyuk may look cumbersome. In constrast, the Yukawa couplings to the 

physical Higgs field h are quite different since the leptons obtain masses from both the bare 

mass term and the Yukawa couplings: 

^Yuk = --mcl (viV3 + WV4) - i-{mA - m3)cese ( ViZV4i? - V^Vsl) 



l^aL ma5afi-{mA-m3)Ul^JjR^p ^jg^ + h.c.j (25) 
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2.5 Constraints on mixing matrices and lepton masses 

For convenience in the next section, we collect here the constraints on Ul^r, nia and m, : 

CI : niT^ce = m^S0 (26) 

C2 : Y^Ul,^UL3a = Y,Ul,^UL4a = (27) 

a a 

C3 : Y.maUli„UMa=0 (28) 
a 

C4 : £m2[/;,„t/^«=0 (29) 
a 

where ? = 1, 2 in C2, C3 and C4. They will be extensively used to improve the apparent 
convergence of the loop integrals and extract the leading terms in the large mass limit of heavy 
leptons. These constraints are exact and can be readily derived. The constraint CI is from 
diagonalization of while C2 represents unitarity of Ul. After rotating the column vector 
y-£ to its third component, the first two columns in the last row of m^. vanish. This yields 
(f/^mjf/2)4i = ('"£)4i = for ?' = 1,2, which is C3. In addition, we find that (m£m^)4/ also 
vanishes for i = I, 2, which gives the last constraint C4. For the sake of notational simplicity, 
we sometimes also use the Latin letters z, j and numbers, which enter through the charged 
current matrices as the indices for the corresponding charged leptons. 



3 Two-loop induced neutrino masses 

Now we calculate the radiative mass of the neutrinos Vi 2 that are massless at tree level. This 
is given by their minus self-energy evaluated at the zero momentum. We thus need to calculate 
the amplitude for the transition, v,x — > vj^, with z, j = 1, 2. There is no contribution at one 
loop. This arises because, while the neutral current does not couple Vi,2 to the massive ones 
V3 4, the charged current involving Vi^2 is purely left-handed and thus cannot induce a mass for 
a massless particle. 

At two loops, we note first that a diagram with at least one of the two external lines con- 
nected to a virtual Z boson cannot contribute. This is because, if it did, removing this virtual Z 
line would also do since Z couples diagonally to Vi^2 and conserves chirality. But this would 
contradict our claim at one loop. The external lines must therefore all connect to virtual 
bosons. Finally, the two external lines cannot connect to the same virtual due to charge 
conservation. This leaves with us the single diagram shown in Fig. 1. 




Figure 1 . Diagram contributing to — zXy,- 
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We shall evaluate the radiative mass in unitarity gauge. We first simplify and classify the 
contributions from the diagram. Then we apply the constraints CI — C4 to reach manifest 
convergence in loop integrals and to get prepared for isolating leading terms in the seesaw 
limit. Finally, the contributions are expressed in terms of some standard parameter integrals. 

To start with, we note that the external Vtj (z, j = 1 , 2) have no right-handed couplings to 
the corresponding virtual charged leptons (-a.fi- The diagram then decomposes into four terms 
according to the chiralities of the two vertices involving the virtual neutrino v^- After some 
algebraic work, we can remove all 7 matrices in favor of the products of loop momenta and 
obtain 



^""J' 4(4;r)4 

where Uij are the spinors for external neutrinos, and ^ gives the radiative neutrino mass. The 
T functions are 

T^^ = m«#^^«^4^M«^M/3i^^''(«,/3;fc) (31) 

where the loop functions F are dimensionless functions of the mass ratios. Upon Wick rotation 
to Euclidian space, they become 



D(a,/3;A:) 

+2(pV -q + lif-q)-- q^P^) -4q^(p ■ q + q^)] (32) 

where the notations are 

D{a,m = [{P + qf + n] [p^ + ra] + 1] [q^ + rp] [q^ + l] 

r.= 4'^« = 4' // = ^/A/A (33) 

Here the summation over the virtual lepton flavors a, /3, is implied in the T functions, and 
is manifestly symmetric as expected for Majorana particles. 

Since only massive neutrinos enter the right-handed charged current, the virtual V]^ in the T 
functions is actually restricted to V3,4. Using the explicit forms of Wl^r shown in eq. (I22l) . the 
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T functions decompose into 

+ V2ceS0{UiAaUuii+UL3aUj^4l3) ^F/^^ + rngF/^^] } 
= 2mampm^^Ul,^UljpUMaUMp[mclF,^''-m3slFf''] 

T^^ = maUl^^UljpUR4a{V2ceseU^3p[Fi'^-Ff'^]+2Ui^p[slF,^^^ 

where for brevity the first two arguments a, /3 of the F functions are suppressed while the 
third one k appears as a subscript 3 or 4. In addition to improving apparent convergence, the 
main merit of applying the constraints CI —C4 is to subtract heavy leptons V4, x from the 
loops. This avoids manifestly in the contributing terms some large numbers that are actually 
balanced by the small matrix elements mixing the light and heavy leptons. Furthermore, this 
facilitates the extraction of the leading terms that can survive upon being multiplied by the 
mixing matrix elements and summing over light flavors a, /3, for which the hierarchical limit 
I ^ ra ^ rj, works very well. We stress that we are not discarding the contributions from 
heavy leptons but are combining them in a judicious manner with those from light leptons 
before numerical analysis is done. In the following subsections we shall reduce the T functions 
using the constraints. 

3.1 Reduction of T^^ 

We note first of all that the numerator of F^^ is separately linear in and q^. Take as 
an example. By decomposing p^ = {p^ + — the first term cancels the corresponding 
factor in Z)(a,/3;fc) so that its contribution to F^^ is independent of a. The constraint C2 then 
implies that it does not survive in T^^ upon summing over a. We can thus effectively set in the 
numerator of F^^, p^ — * — and similarly ^ —r^: 

F^(a,P:k) -|4+.„.,-4(.„ + .,)]//^^^, (35) 

where the arrow means equality when multiplied by U factors and summing over a, /3. To 
go further, we have to cope separately with the four terms in T^^ according to the Ui factors 
involved: 

+ V2 [UiAaUL3fi + UL^aUtAH Tit] } (36) 

with obvious definitions on T^^ etc by comparing with eq. (|34|) . 

Although the first term, T^^, is already convergent upon applying CI due to the subtraction 
between and Fj'^, we can do better by subtracting explicitly the contribution from the heavy 
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charged lepton X- The trick is that, for a term in F^^ that is not proportional to ra we make the 
substitution 

1 1 1 



da{p) 



while for a term that is proportional to ra, we do as follows 

fa fa fA 



P^ + fa 



P +fa P + fA 



ea{p) 



(37) 



(38) 



The legitimacy of the substitutions is guaranteed by the constraint C2. Thus, 



'T'LL 
^33 



p-q 



rdsip + q) 



[p2 + l][^2+i] 

X [4da{p)dii{q) + ea{p)ep{q)-4eaip)dp{q)-4da{p)eii{q)] (39) 

The second term, T^^, is multiplied by UiAaUiAp so that we have a choice of whether to 
use the constraint C2 (i.e., eq. (l37l) ) or C4 (eq. (l38l) ) for the terms proportional to or r^. It 
turns out that the latter is better as it can reduce the amount of work by bringing down more 
factors of r^^^ for light leptons a, /3, which makes the corresponding term subdominant in the 
hierarchical limit. The last two terms may be similarly manipulated. The results are summarized 
as follows: 

p-q 



I44 



[4 + r«r^-4(r« + r^)] JJ ^-^-l:A__daip)dpiq) 



X 



m4c| 



rLL 
'43 



_{p + q)^ + f3 {p + qY + fA_ 
[p2+^][^2+l] - 1)^/3 (^) + {^-fa)ef,{q)\ da{p) 



X 



m4 



+ ■ 



m3 



cese 



(40) 



_{p + q)^ + fA {p + qY + f3. 

/3 and i ^ j. Since a, /3 are summed over, this amounts 



while T^^^ is obtained from T^^ by a 
to symmetrizing T^^ in z, j. 

The advantage of the above results can be understood by recalling that we now only need to 
sum over light flavors a, /3 in T^^. Since 1 ^ ra^p ^ r3, it is numerically very good to set = 
fji — f3 = 0. For instance, the largest r-r ~ 5 x 10^^ while r3 ~ 6 x 10^24 fQj- ~ 0.2 eV. This 
will not introduce mass singularities in the loop integrals. In addition, when a term proportional 
to m3 is accompanied by one proportional to m^, we ignore the former since it cannot make a 
significant contribution to the radiative mass. (Note that T^^ is exceptional since m^c^ = m^s^Q.) 
Although the above argument is self-evident, we have inspected and compared carefully all of 
the terms to verify it. This simplifies considerably the integrals to compute: 

• m34r2{4[^2(0)-^2(r4)] + 8[=ri(0)-=ri(r4)] + ^o} 
-m/^CQArl^^2{fA) 

,43 ^ -m4Ce5e4r|[^2(r4) + ^l(r4)] (41) 

and T-^l" = T^^, where the loop integrals ^ are defined in Appendix A. These functions are 
independent of a, /3 and depend only on r4. 



T 

^33 

7-'L>L> 
L 

Z 



44 

LL 
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3.2 Reduction of T 



The second term in T^^ is doubly suppressed by ^35^ compared to the first one and will be 
ignored from the start. Since the numerator in the integrand of is again linear in and 
q^, they may be replaced by — r« and —r^ respectively employing the constraint C3. For the 
2{p-qY term in the numerator, we decompose as follows, 

2{p .qf=p.q{[{p + qf + r^] - [p^ + r«] - [q^ + rp] + [r« + - r4]) 

The first term is cancelled by the same factor in the denominator D making the integrand odd 
in p, and thus vanishes upon integration. The second term again cancels a same factor from D 
and is killed upon summing over a by the constraint C3, and the same happens with the third 
term as well. The numerator now becomes effectively, 

-{S + r4p ■ q) + {p ■ q + 2){ra + rp) - rarp 

Now we make the substitutions in eqs. (l37l38l) as we did in the previous subsection, though 
employing this time the constraint C3, to obtain, 

+{p-q + 2) [ea{p)dji{q) + ep{p)da{q)] -ea{p)ep{q)} (42) 
Since it is now legitimate to sum only over light flavors a, j8, the above simplifies to 

if ^ -rl{S'3^2ir4) +4^^i(r4) +^0(^4) +^4^2(^4) +2jri(r4)} (43) 
where the new integrals '3^ are also defined in Appendix A. 

3.3 Reduction of T^^ 

This chirality-mixed part from the two vertices involving the virtual neutrino contains the 
most number of terms in F^^: 

= maUl,^UljpUMa{V2ceseUt^3p [Fi'^-Fi'^] +2clUt^pFi'^] (44) 

where we have dropped the s^Ff^ term as one cannot rely on it to induce a reasonable mass 
due to a tiny ~ 10^^^ at ~ 0.2 eV and 1714 ~ 200 GeV, for instance. 

The numerator of the integrand in F^^ is linear in p^, which can thus be replaced by — r« 
using the constrain C3. On the other hand, since the numerator is quadratic in q^, we must 
distinguish between the two terms in T^^ which are proportional to Ui2j5 and Ui4p respectively. 
For the first one, we can only set one factor of q^ to —rp using C2. After this, we apply C2 and 
C3 via the substitutions in eqs. (I37I38I) and obtain, 

Fi^-F,^^ jjj^^^^^^[{p-q + q^ + 2)ea{p)ep{q) + 2p.qea{p)dp{q) 

-A{p-q + q^+l)da (p) ep{q) +4 {p ■ q - {p ■ q) ^) da {p)dp (q)] (45) 
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The summation over light flavors a, j8 then yields the result in terms of the standard integrals: 
F^^-Ff"^ ^ r||r4[^o + 4^1] + ^0 + 6[ JTi (0) - ^i(r4)] 

+A[^2{0) - ^2(/-4)] - 2r4 Jr2(r4) + [%{0) - } (46) 

For the second term proportional to Ujj,^ in T^^, we can set two factors of to —r^ because 
of C2 and C4. The subsequent manipulation based on the constraints and eqs. (|37I38I) is similar, 
and gives 

Ff^ ^ r2{r^(^ro(r4) + [jri(r4)+2^^i(r4)]+4^ri(r4)+4[^^2(r4) + ^2(r4)]) 

-4^2(^4) -2[^l(r4)+r4Jr2(r4)]} (47) 

where the terms suppressed by will be ignored from now on. 

To finish this section, we summarize the terms in the radiative neutrino mass as follows: 

+2^rarpUR4aUR4pm4clFf^ 

+V2C0S0 {maUR4aUL3p+mpUR4iiUL3a) W^-if^) 

+24 {maUR4aUi^p + mp Ur4p Uiao) } (48) 

where the relevant functions are given in eqs. (HTl l43l l46l |47|) in terms of the standard inte- 
grals calculated in Appendix A. The summation over the light charged leptons ia and ip is 
understood in the above. 



4 Numerical analysis 

Now we investigate whether we can accommodate the neutrino masses measured in oscillation 
experiments. Our starting formula was given in (|48l) which involves the light-heavy mixing 
parameters in addition to the upper-left 3x3 submatrix of Ut. From eq. (l22l) we see that 
the latter is just the leptonic mixing matrix measured in oscillation experiments to very good 
precision. However it is no more exactly unitary, and the deviation from unitarity is determined 
by the light-heavy mixing. A realistic numerical estimate should take all this into account to 
avoid a misleading conclusion. Although a global fitting to the lepton mixing parameters is 
possible with radiative corrections included, our main result on the seesaw scale required to 
reproduce the neutrino masses is independent of this fitting. 

Both matrices m^m^ and m^m^ for the charged leptons have the hierarchical structure 

M^(Jtl), (49, 
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where B and d are respectively a 3 x 3 and 3x1 matrix, whose entries are much smaller in 
magnitude than the positive number A. Then, the submatrix of the diagonalization matrix that 
mixes the small and large entries can be estimated as K ^ dA^^. Application of this to niEin'^ 
and m\mE yields for a — jl, r: 

U^a ~ (msM) = (r3/r4) 1/^ = 

UR4a ~ {m3/m4y/^{ma/m4) = e{ra/r4y/^ (50) 
where (fT8l) is used. And the unitarity violation in the submatrix of light leptons is, for i = 1 , 2, 

L Ul,aUL3a = -Uli^UL3x-e^ (51) 
a=e,fi,r 

Consider first the case in which m4 is not very large. This is the range of parameters that 
is particularly relevant to LHC physics. A heavy active lepton, especially the charged one 
is supposed to be accessible if it is not much heavier than several hundred GeV. Our estimate 
of heavy-light mixing parameters is still good enough since m4 is much larger than the light 
lepton masses. Using the estimates in eqs.(fT8l [50l) (but not yet the one in (|5TI)). we find that 
the three classes of contributions to in eq. (l48l) consist of the following terms in units of 

RR : rarpUl^^Ul^p 

RL : ra{Ul,^Uljp + Ulj^Ul^p)UE3p,raiUl,^Ul^^ (52) 

where each term is to be multiplied by a coefficient that is a sum of integrals as can be obtained 
from eqs. (147114311461 l47l) . The point is that these coefficients are order one numbers for r4 not 
very large. Then, independently of the mixing matrix of light leptons, it is safe to say that 

\^ji\<l.SxlO-^m3 (53) 

Since no light neutrinos can be heavier than an eV from cosmological considerations, there is 
no hope to induce a large enough radiative mass mi or from m^. Therefore, the minimal type 
III seesaw model cannot accommodate oscillation data if the heavy leptons have an intermediate 
mass. To put another way, the oscillation data already excludes the possibility that the active 
heavy leptons in the model would be accessible at LHC. 

It is interesting to ask whether there is a chance at all in the model to induce a large enough 
neutrino mass. For this purpose, we study the seesaw limit in which 1114 blows up. Then is 
a sum of the following terms (again in units of 2^^ n^^rn^Gj^m^): 

LL : r3r4[^o]-Ul^^UljpS[rl^2{r4)]-{Uha + Ulja)^ 
RR : -raUli^rpUljp2[r4%{r4)+rl^2{r4)+2r4^l{r4)] 

RL : [raUli^y^4[r4% + ^o]-raUl,^Uljp4[r4^,{r4)+rl^^^ 

+ (54) 



13 



All combinations of loop integrals in the square brackets are 0(1) constants up to logarithmic 
corrections in the large r4 limit. We have also taken into account the unitarity violation estimated 
in eq. (ISTI) . Because of the estimates employed, the relative sign and factors of two between 
terms in the above cannot be taken seriously. But this does not preclude us from making a 
definite conclusion as shown below. 

To induce a mass of 0{m'i), some terms in eq. (|54l) must be above 10^. This obviously 
requires a large r^. But even this is insufficient. On the one hand, the terms not multiplied by 
r4 factors outside the square brackets can be safely ignored; on the other, all remaining terms 
are controlled by r^r^. We must therefore require r-^r^ ^ 1. This corresponds to the combined 
limit in terms of the original parameters in Lagrangian, ^ vr^ ^ mw- In the limit, only the 
first term in the LL class is relevant: 

J^j, ~ 2-^K-Sn%Glm^r^r^ [^o] (55) 

Inspection of our derivation shows that this is the term that is doubly suppressed by unitarity 
violation between the third row and the first two rows of the light lepton mixing matrix. But 
unfortunately it is impractical to measure the violation down to the level that we are interested 
in, i.e., ~ 0^ = r3/r4 = m^/m\. The information on the indices (z, j) is lost also because of 
the estimates employed. This means in passing that our analysis on the neutrino masses in the 
above limit is independent of a detailed fitting to the leptonic mixing parameters. We find it 
is natural for the model to favor the normal hierarchy scenario; namely, a larger m-i seeds a 
smaller mi 2- For the purpose of illustration, we assume m\ = 0. The solar and atmospherical 
oscillation data then give m2 ~ 8.7 X 10^3 ev and m3 ^ 4.9 X 10 2 eV respectively, which can 
be fulfilled by requiring 

m4 ~ 4 X 10^^ GeV (56) 
This is roughly the scale of grand unification. 

5 Conclusion 

The minimal type III seesaw model introduces a lepton triplet on top of the particles in SM. Two 
neutrinos out of four are massless at the tree level, but they are not protected by any symmetry 
from getting a radiative mass at the quantum level. We have shown that the latter takes place 
first at two loops, and determined it in terms of some parameter functions. By employing 
realistic estimates of the mixing parameters between the light and heavy leptons, we studied the 
pattern of the neutrino masses. We found that it is not possible to accommodate the spectrum 
determined in oscillation experiments if the heavy leptons have a mass that would be within the 
reach of LHC. However, if the seesaw scale is as large as that of grand unification, it is still 
possible to accommodate the spectrum in a nice manner: one light neutrino gets mass directly 
from seesaw while the other two get a radiative mass. The model would then contain nothing 
new but the tiny neutrino masses. The main message extracted from this work is therefore, if 



14 



LHC sees something like a triplet lepton, it definitely comes from a structure that goes beyond 
the economical one as originally suggested. 
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Appendix A: Loop integrals 

The loop integrals in the final result of T^^ (see eg. ((41)) ) are defined as 



^0 



p-q 



Di{r)p^q^ 

p-q 
p-qrA 

Di{r4){p + q)^ 



where 



Di (r) = [p2 + ,4] [p2 ^ 1] [^2 ^ [^2^1] + qf + r] 
The new integrals appearing in T^^ and T^^ are respectively, 



1 



1 

1 



and 



II D2 



^1 



Djq' 



where 



(57) 



(58) 



(59) 



D2 = {p + qf[{p + q)^ + r4] [q^+l] [q^ + [p^ + r4] 
There is another integral in calculating T^^ that can be related to those already defined: 



(60) 



(61) 



II Di(r 



The basic technique to compute the above integrals is to use fractions and the one-loop 
integrals in « = 4 — 2e dimensions: 



{Anf 



d"p 



1 



{2nY [{p + q)^ + r][p^ + a\ 
d"p p ■ q 

{2k)" [{p + q)^ + r][p^ + a] 



= {Anfq^ 



r(e) — I dx \ng{a,r) 
Jo 



1 /■! 
— -r(e)+ / dx x\ng{a,r) 



(62) 
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where 



g{a,r) = q^x{l —x) +rx + a{l —x) (63) 
Introducing the abbreviations, 

g(a,r) =x{l -x){l -y) + [rx + a{l -x)]y 

gQ = —( V' <?H^^ = ~( V 

^(r4,r4) gyr^.r) 

^^^^^lng(r4r)^lng(l,r)_^lng(0,r) ^^^^ 

(r4-l)r4 l-r4 r4 

and denoting the parameter integrals in the form, 

X= C dx [^dyI[X] (65) 
Jo Jo 

where X enumerates all of the defined integrals, the integrands are 

jc(l-y)^(r) 



y{\ -y + r4y) 



r4- 1 y{l -y + r4y) 

^H^o = 7^7^ rln^^^ (66) 

r4-ly^{l-y + r4y) gi{r4) 



for the sequence, and 



1 -y + r4y 



imir)] = j^-. ^lnli(r) 

1 -r4 1 -y + r4y 

imr)] = rr-^l"^i('-) 

l-r4y{l -y + r4y) 

im = -- J ngo 
r4y{l -y + r4y) 

im = Inio (67) 

r4 1 -y + r4y 

for the ^ and ^ sequences. 

The above integrals have a magnitude of order one or smaller for r4 not very large, and can 
be readily integrated numerically. This is a sufficient message for the first part of our numerical 
analysis in section IH For the analysis in the heavy mass limit, we need the leading terms of the 
integrals. We obtain them in two ways. One is to use the techniques and formulae developed 
already in the literature ll34l [30l . and extend them slightly to cover all cases occurring in our 
integrals. (There is a typographic error in expansion (ii) on page 230 in Ref [[34ll : ^\y? a should 
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have a plus sign instead of minus.) The leading terms can also be extracted directly. For 
illustration, we calculate below the integrals ^\{ra) and .^(r4) that appear most frequently in 
eq. (l54l) . We finish first the integration over 3; in terms of logarithm and dilogarithm functions 
using 

Kb) = ^—H^ + {b-l)y] = -U2{l-b) 
Jo y 

^1 \n[l + {b-l)y] 



J{b,r) = (r-1) / dy- 
Jo 



l + {r-l)y 



= uJ^^-\ VLi^f^Vln^lnr + lln^r (68) 
\{b-\)r) \b-l) b-l 2 ^ ' 

where > 1 , r > 1 . Denoting 

^ r^x+\-x ^ ^4 ^ ^4 ^ ^ bj - r4 ^^^^ 

x{\—x)' x(l— jc)' 1—x' bi—\ 

with Z72 > ^1 > ^3 > ''4 > 1 > > for X G (0, 1), and using the abbreviations 

I, = I{b,),J,=J{b,,r4) (70) 
we express the integrals as follows: 

(r4-l)jri(r4) = lxdx[{[h-l2]-[Ji-J2])-^j^} 

r4{r4 - l)Mr4) = J xdx^r^^h- h] - [Ji- Jl]) - {r4- I) {[h- h] - [J2-J3]) 

J2—J1 J2 — J3 



(71) 



r4 - 1 r4 

Since none of and 7, diverges as a power as becomes large, we have for r4 ^ 1 , 

(r4-l)^i(r4) = j xdx[[h-l2]-[Ji-J2]}+0{r^') 

r^^r^ - 1)^2{h) = I xdx^r^iik-h] - [J3-J1]) + [h-h] + [^1 -^2] } + 0{r^'){12) 

To extract the leading terms, we have to expand the first combination in .^2(^4) to 0{r^^) 
and all others to 0(1). Consider the latter first. Since all bi ^ 1 for r4 ^ 1, we use Landen 
identity of dilogarithm for the last two combinations in ^2' 

h-h = ^ln(M3)ln^+Li2(l-^2')-Li2(l-^^^) 

= [21nr4-lnx-21n(l - x)]\nx + 0{r^^) 

Ji-h = [Li2(ai/r4) -Li2(ai) +ln(Z?i - l)lnr4] - (ai ^ a2;^i -> ^2) 

= Li2(l-;c(l-;c))-Li2(;c)+lnxlnr4 + C>(r4^) (73) 
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Then 



Bi = j xdx{[h-h] + [Ji-J2\} 



1 lllt^ 1 



+ 4v^i(l/6) + iv^i(l/3) + 0(r4i; 



2 36 ' ' \2 

1 



f« 0.435 + 0(r4^) (74) 

where V^i(z) = Inr(z) is the trigamma function. Since 

[h -h] - [J I -Ji] = -[h-h] - [J I -J2]+0{r^'), (75) 
this also gives the leading term 

r4^i(r4) = -52 + 0(^4 (76) 

The first combination in ^ is more complicated. Using Landen identity and expansions of 
Li2(z) at z = and z = 1^, we have 

r4 X I —X 

1 1 —X 

J3-J1 = Li2(ai)- 02(03) lnr4 + 0{r^^) (77) 

^4 X 

The combination is thus 

= nj xdx{[h-h]-[J3-h]} 

= J dx{{\-x)\n{\-x)+r4x[U2{a3)-U2{ai)]} + 0{r-^) (78) 

The Li2(a,) terms can be worked out by integration by parts, noting that ai ,3 = 1 at x = 1 while 
ai = 1 and as = at jc = 0: 

J xdxU2{ai^3) = ^ + ^ J 6fxx^ln(l -ai,3) ^^^^^'^ (79) 

^ JLi2(z) In(l-z) . ^ ^- u • ^ • -1 

where = is applied. Upon expanding the integrand in r . , we arrive at 

dz z 

J X(ijc[Li2(a3) — Li2(ai)] = - — J dx 



X 



4^-1^ + ^^^^") '''' 



so that 5i = 1 - — + 0{r^ 
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We collect below the leading terms for all integrals. 



and 



^0 = ^-^Co + 0(r4^) 
1 



r4^i(0) = l-^ + C>(r4i) 

r4^i(r4) = ^ + ^-^Co + 0(r4i) 

r|^2(0) = -l + ^;r2-lnr4 + 0(r4i) 

rl^2(r4) = \-'^^\c^ + 0{r-') (81) 



;r2 



r4^o(r4) = -y+Co + 0(r4i) 

r2^i(0) = l-^+lnr4 + ^ln2r4 + 0(r4-i) 
rl^x{n) = 3-Co + \nr4 + 0{r^^) 

rim) = y + 0(r4i) 

rl^A) = -7 + y + Co-21nr4 + ln2r4 + 0(r4i) (82) 

D 

= 3-Co + lnr4 + 0{r^^) (83) 

with Co = 2V3Cl(;r/3) = -4;rV9 + l/6v/i(l/6) + l/6v/i(l/3) ^ 3.51586, where CI is the 
Clausen function. These leading terms have been numerically verified. 
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